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1.  Introduction 

Nitrate  esters,  such  as  nitrocellulose  (1)  and  nitroglycerin  (2),  are  one  of  the  major  classes 
of  energetic  compounds,  along  with  nitramines,  nitroaromatics,  nitroheterocycles,  etc.  [1,2]. 
Nitrocellulose  and  nitroglycerin,  for  example,  are  frequently  components  of  explosive  and 
propellant  formulations.  Accordingly,  any  effort  to  quantitatively  analyze  and  model  the 
combustion  behavior  of  such  formulations  requires  a  knowledge  of  nitrate  ester  decomposition 
processes  and  the  associated  energetics.  As  an  initial  step  toward  acquiring  such  data,  we  have 
carried  out  a  computational  study  of  several  decomposition  routes  for  the  simplest  nitrate  ester, 
methyl  nitrate  (3). 
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2.  Method 

We  have  used  two  density  functional  procedures,  as  incorporated  in  the  program  Gaussian 
94  [3]:  (a)  The  exchange  and  correlation  functionals  were  the  Becke-3  (B3)  and  the  Lee,  Yang  and 
Parr  (LYP),  respectively  [4,5],  and  the  basis  set  was  the  6-31+G(d).  (b)  The  LYP  functional  was 
replaced  by  the  Perdew-Wang  (PW91)  [6],  and  the  6-31G(d,p)  basis  set  was  used. 

3.  Results 


(3): 


We  have  investigated  three  possible  decomposition  pathways,  shown  as  eqs.  (1),  (2)  and 
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In  eqs.  (IB)  and  (3),  TS1  and  TS2  are  transition  states. 

The  optimized  geometries  of  the  reactants,  products  and  transition  states  in  eqs.  (1)  -  (3)  are 
given  in  Table  1.  The  two  sets  of  calculated  structures  are  overall  in  very  good  agreement  with 
each  other  and  with  the  available  experimental  data.  In  Table  2  are  the  total  and  zero-point 
energies.  We  found  that  the  reaction  in  eq.  (3)  requires  an  initial  rotation  around  the  C-0  bond,  as 
shown  in  eq.  (4): 
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The  geometries  and  energies  of  both  the  ground  state  3A  and  its  conformer  3B  are  given  in  Tables 
1  and  2. 

In  Figure  1  are  presented  the  relative  energies  at  0  K,  including  zero-point  contributions,  of 
the  various  stages  in  the  processes  shown  in  eqs.  (1),  (2)  and  (3).  The  corresponding  enthalpies  at 
298  K  are  given  in  parentheses  [7].  In  Table  3  are  compared  the  calculated  and  experimental 
enthalpies  for  the  two  dissociation  processes,  eqs.  (1A)  and  (2),  and  the  overall  AH  for  the 
reactions  in  eqs.  (1)  and  (3). 

4.  Discussion 

The  data  in  Table  3  present  an  interesting  contrast.  The  overall  AH  for  the  reactions  going 
through  transition  states,  eqs.  (1)  and  (3),  is  given  reasonably  accurately  by  the  B3/LYP 
procedure,  better  than  by  the  B3/PW91;  however  the  reverse  is  true  for  the  dissociation  processes, 
eqs.  (1A)  and  (2),  for  which  the  B3/LYP  values  are  too  low  by  about  7  kcal/mole.  It  has  already 
been  noted  in  the  past  that  the  B3/LYP  combination  tends  to  underestimate  the  energy  requirements 


for  the  detachment  of  nitro  groups  [8,9].  This  problem  is  not  overcome  by  using  a  larger  basis  set, 
as  can  be  seen  in  Mebel  et  al  [8]  and  as  we  now  confirmed  by  repeating  the  B3/LYP  calculations 
for  eqs.  (1A)  and  (2)  at  the  6-31 1+G(d,p)  and  6-3 1 1+G(2df,p)  levels  with  no  significant 
improvement.  However  we  have  shown  earlier  [9],  for  the  specific  cases  of  H3C-NO2  and 
H2N-NO2,  that  better  results  can  be  obtained  with  the  B3/PW91  combination,  and  our  present 
experience  reinforces  that  conclusion,  for  splitting  off  ONO2  as  well  as  NO2. 

Figure  1  shows  that  the  two  computational  methods  are  in  good  agreement  concerning  the 
activation  barriers  to  forming  transition  states  TS1  and  TS2.  The  reliability  of  B3/LYP  activation 
energies,  in  particular,  has  been  investigated  extensively  [10-16];  while  an  occasional  problem  is 
encountered  [11],  in  general  the  results  are  quite  satisfactory. 

With  regard  to  the  three  methyl  nitrate  decomposition  routes  that  have  been  investigated,  we 
conclude  therefore  that  the  energy  requirement  for  the  first  step  in  eq.  (1)  is  comparable  to  that  for 
eq.  (3).  Accordingly  both  of  these  processes  can  be  expected  to  occur,  although  the  second  step  in 
eq.  (1),  leading  to  TS1,  has  an  additional  significant  energy  barrier.  It  should  be  noted  that  eq.  (3) 
is  the  only  one  of  the  three  pathways  that  involves  a  net  release  of  energy. 

5.  Conclusions 

We  find  two  likely  initial  decomposition  paths  for  methyl  nitrate  to  be  (a)  loss  of  NO2, 
followed  by  eventual  rearrangement  to  H2C-OH,  and  (b)  formation  of  H2C=0  and  HONO.  The 
initial  energy  requirement  for  each  process  is  about  40  kcal/mole,  and  the  second  is  exothermic, 
with  AH  (298  K,  experimental)  =  -16  kcal/mole  of  H3C-ONO2.  From  the  standpoint  of 
computational  methodology,  both  the  B3/LYP  and  the  B3/PW91  density  functional  procedures 
were  overall  effective;  however  it  was  again  observed  that  the  B3/LYP  tends  to  underestimate 
dissociation  energies  involving  the  loss  of  NO2  (and  ONO2  as  well). 
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Table  1.  Optimized  geometries.3’13 
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1.426,  1.411 
1.177,  1.177 

(0.958) 

(1.432) 

(1.170) 

H-Oa-N:  103.0,  102.5 
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H— C-Oa:  96.6,95.8 

aFor  each  distance  or  angle,  the  first  entry  is  the  B3/LYP  result,  the  second  is  the  B3/PW91. 
Experimental  data  are  in  parentheses.  They  are  taken  from  M.  D.  Harmony,  V.  W.  Laurie,  R.  L. 
Kuczkowski,  R.  H.  Schwendeman,  D.  A.  Ramsay,  F.  J.  Lovas,  W.  J.  Lafferty  and  A.  G.  Maki, 
J.  Phys.  Chem.  Ref.  Data  8,  619  (1979). 
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Table  2.  Calculated  total  and  zero-point  energies. 


Molecule 

Total  energy,  hartrees 

Zero-point  energy,  kcal/mole 
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-114.50884 

-114.45657 

16.8 

16.8 

H2C-OH 

-115.06294 

-115.01734 

23.4 
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18.7 

HONO 

-205.70966 

-205.62186 
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Table  3.  Comparison  of  calculated  and  experimental  AH  values,  at  298  K,  for  reactions  shown  in 


eqs.  (1A),  (1),  (2)  and  (3). 


Reaction 

AH,  298  K  (kcal/mole  of  H3C-ONO2) 
B3/LYP  B3/PW91  Experimental3 

1A:  H3C-0N02  — 

h3c-o 

+  no2 

33.9 

38.4 

36.5, b  41, c  40.7d 

1:  H3C-0N02  — 

->-h2c-oh 

+  no2 

32.2 

34.3 

30.3d 

2:  H3C-0N02  — 

H3C  + 

no3 

74.3 

76.4 

80.9, b  80.7d 

3:  H3C-0N02  — 

h2c=o 

+  HONO 

-12.3 

-6.8 

-15.7, b  -16c 

aObtained  from  experimental  heats  of  formation. 
bJ.  J.  P.  Stewart,  J.  Comp.  Chem.  10,  221  (1989). 

CS.  G.  Lias,  J.  E.  Bartmess,  J.  F.  Liebman,  J.  L.  Holmes,  R.  D.  Levin  and  W.  G.  Mallard,  J. 
Phys.  Chem.  Ref.  Data  17,  suppl.  1  (1988). 

dD.  F.  McMillen  and  D.  M.  Golden,  Ann.  Rev.  Phys.  Chem.  33,  493  (1982). 
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Figure  1 .  Relative  energies  at  0  K  and  relative  enthalpies  at  298  K  (in  parentheses)  of  the  various 
stages  in  the  decomposition  processes  shown  in  eqs.  (1),  (2)  and  (3).  The  entries  at  the 
left  are  B3/LYP  results;  those  at  the  right  are  B3/PW91.  All  data  are  in  kcal/mole  of 
H3C-0N02. 


